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Abstract. 

X-ray reflectometry (XRR) provides researchers and manufacturers with a non- 
destructive way to determine thickness, roughness, and density of thin films 
deposited on smooth substrates. Due to the nested nature of equations modeling 
this phenomenon, the inter-relation between instrument alignment and parameter 
estimation accuracy is somewhat opaque. In this study, wc intentionally shift incident 
angle information contained in a high-quality XRR data set and refine this shifted 
data using an identical structural model to assess the effect angle misalignment has 
on parameter estimation. We develop a series of calibration curves relating data 
misalignment to variation with layer thickness and density for a multilayer GaAs/AlAs 
Certified Reference Material on a GaAs substrate. We then test the validity and 
robustness of several approaches of using known thickness and density parameters 
from this structure to calibrate instrument alignment. We find the highest sensitivity 
to and linearity with measurement misalignment from buried AlAs and GaAs layers, 
in contrast to the surface layers, which show the most instability. This is a fortuitous 
result, as buried AlAs and GaAs exhibit the highest long-term stability in thickness. 
Therefore, it is indeed possible to use reference thickness estimates to validate XRR 
angle alignment accuracy. Buried layer mass density information also shows promise 
as a robust calibration approach. This is suprising, as electron density of buried layers 
is both a highly-correlated phenomenon, and a subtle component in the XRR model. 
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31 1. Introduction 

32 X-ray Reflectivity (XRR) has been used to measure thin films since its discovery by 

33 Kiessig in 1931 [1]. Parratt in 1954 [2] developed a first-principles approach to modeling 

34 XRR data. Since that time, XRR has matured as the preferred method for non- 
35 destructive evaluation of thin film thicknesses (e.g., Lekner, [3] Chason and Dalliant 

36 j5]). XRR is an interference phenomenon between layers of distinct electron density. 

37 This interference phenomenon manifests as fringes in reflection intensity as a function 

38 of incident angle, 9. Its ability to detect interfaces non-destructively is simultaneously 

39 its strength and its weakness as a technique: layers with diffuse interfaces, such as 

40 inter-diffusion layers and layers with high roughness, will often be difficult, if not 

41 impossible to model. Because reflection takes place at glancing incidence, the technique 

42 also typically illuminates a large sample area during measurement, it is difficult to 

43 analyze the thickness of films with large-scale lateral variation. However, in the realm 

44 of smooth films of high lateral uniformity (such as optics coatings and semiconductor 

45 manufacture) XRR has become an invaluable tool. One significant advantage to XRR 

46 over visible-light thickness determination methods (such as optical refiectometry and 

47 ellipsometry [6]) is that the index of refraction, in the hard x-ray spectrum, is near 

48 unity for all materials. For visible light thickness determination methods, materials 

49 exhibit large changes in index of refraction, which cause first order changes to thickness 

50 determination. For XRR, these corrections due to index are negligible, making XRR 

51 an ideal approach for International System of Units (Sl)-traceable measurements of 

52 thickness. For accuracy in thickness determination, we need calibration artifacts for 

53 inter-comparison measurements across XRR tools manufactured and used across the 

54 world. 

55 In the last decade, an international collaboration under the Versailles Project on 

56 Advanced Materials and Standards (VAMAS) has performed round-robin studies on 

57 several candidate materials to be used as thickness standards for XRR [71 [8]. One 

58 structure studied by VAMAS, a three-repeat bi-layer of GaAs/AlAs (total of six layers) 

59 deposited on a GaAs wafer, was developed by the National Metrology Institute of Japan 

60 (NMIJ) as a pre-standard for a NMIJ Certified Reference Material (CRM). Data taken 

61 in 2004 by researchers at NMIJ on a similarly deposited structure, is the focus of this 

62 study. The National Institute of Standards and Technology (NIST) has been developing 

63 Bayesian approaches (see Sivia f5]) to estimate uncertainty in modeling XRR parameters 

64 [MM- 

65 2. Experimental Details 

66 XRR data are collected as a series of incident angle, 9\ and refiected intensity, I^, data 

67 pairs, {61, iP) stepped over points in a range, / = 1, 2, . . . , from a starting incident 

68 angle, ^start ending incident angle 6**^^. These data should represent the specular 

69 refiection from the surface of the sample. Specular reflection occurs when the incident 
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70 and reflected angles are equal: 6"^ = 6^'. Our model for XRR starts with the assumption 

71 of specular reflection. However, instrumentation effects, such as sample misalignment, 

72 6, may introduce inequality to this equation. An XRR measurment instrument does not 

73 directly measure 6^ or 6^ . Instead, a parallel beam of X-rays impinges on a sample which 

74 is rotated to a sample angle, u. Intensity data is collected from a detector which is in 

75 turn rotated to a a detector angle, 9^. We define an aligned instrument (and sample) 

76 to be one, where, u = implies that O'^ = 0. Specular condition is obtained only when 

77 Lo = 6' and 6'^ = 2 X Lo = 6^ + 6''] also called a 6 - 26 scan. 

78 The GaAs/AlAs multilayer structure used in this study was sample 1-08 of the 

79 BAAA4002c series provided by NMIJ (referred to here as 1-08). 



3. Theory 

3.1. Reflectometry modeling 

The first-principles treatment of XRR modeling is discussed in exhaustive detail 
elsewhere (see [3j). However, a brief discussion of several key features of the model 
are useful to understand misalignment impact. XRR measurement modeling typically 
uses a Fresnel homogeneous layer model to treat thin films (j = 1...N) in a stack as slabs 
of a fixed, complex refractive index, nj, with a real component, 6j, and an imaginary 
component, Pj, as defined by nj = 1 — 6j + i(3j. Parratt developed this analysis for 
XRR by solving Snell's Law of refraction for each successive layer in a stack. [2] Reflected 
intensity from any layer can be calculated by starting with the perpendicular component 
of the wave vector from a material layer j of fixed n: 



k,j = k^n] - cos"^ 9^ (1) 

At each interface in a stack of such fixed index, continuous layers, between layers 
j-to-j + 1, we have the following, Fresnel relationship: 

where wave vector, k = 2tt/X, and A is the X-ray wavelength. 

Parratt 's developed a recursion relation for the reflection coefficient for every 
interface, j to j + 1 in a stack, 

1 + ^ ij+i^i+i 
where zj is the thickness of layer j . 

All successive reflection coefficients from Xj_i to Xq can be solved for using 
successive substitutions of Eqns. [2]and|3j with the special case that, since the substrate 
is assumed to be infinite, we have no reflection from its lower interface: X^r jv+i = for 
an N-layer stack. The measured reflection intensity is then = |Xo|^. 

Layer thickness information is "locked" into the recursion though each instance of 
Zj, and layer density information enters as a function derived from the index of refraction. 
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106 rij, for each layer within the recursion. An XRR structural model is defined by a set of 

107 paramters for layer thickness tj = Zj, layer density, pj = f{nj), and layer roughness, aj. 

108 The fundamental observation for this work is that is a function of ^* through each 

109 instance of kzj in the recursive equation. Separation of these parameters is impossible 
no due to the nested, interdependent nature. The focus of this study will be to determine 

111 both severity and parameter refinement trends given intentional shifting of the observed 

112 9^ for high quality data and a model test system. 

113 3.2. Data refinement approach 

114 Our data refinement approach follows the structure and notation of Wormington [12] . 

115 A highly constrained model for calculating simulated measurement pairs, (^;, J/^"^'^^'^), 

116 is then fit to the data and an error between data and calculated intensities is calculated. 

117 It is common to show ot goodness-of-fit (GOF) plots as a function of parameter 

118 variation to establish best fit between data and mathematical model, for example 

119 see [13]. One of Wormington's early achievements with XRR data analysis was the 

120 implementation of alternative strategies for fitting to data. XRR measurements often 

121 scale over many orders of magnitude requiring modification to traditional data weighting 

122 schemes. Wormington implemented several cost functions for XRR data refinement. The 

123 Mean Square Error of the log of refiected intensity, log Ir, values provides a robust cost 

124 function for refining XRR data: 

MSEi,,{p) = Yl [^"^alf - loglf-'^'^iOf, p)] ' (4) 

i=i 

126 MSEiog allows for oscillation data contained several decades below the incident 

127 intensity,Jo, to still contribute to the model refinement. Figure[l]a) shows specular XRR 

128 data, with log of refiected intensity, log/^, in counts as a function of incident /refiection 

129 angle 6 in radians. Note that oscillation data are present over nearly seven orders of 

130 magnitude in intensity. Traditional would ignore contributions form the last five 

131 orders of magnitude of the data. 

132 Wormington's other revolutionary contribution was the use of Differential Evolution 

133 (DE), or more frequently referred to as a Genetic Algorithm (GA), to refine XRR data. 

134 The GA method allows for a large population of potential parameter solutions, p, to be 

135 simultaneously tested against XRR data with one of the modified cost functions, (e.g., 

136 MSEiog). Of this large population of solutions, only the top half "fittest" or highest 

137 goodness-of-fit solutions are allowed to continue in the analysis process. These solutions 

138 are then bred or allowed to intermix parameters p{xj), and then produce a new large 

139 population of solutions for testing against the XRR data. This culling, breeding, culling 

140 model is cycled through hundreds to thousands of populations to generate the global 

141 minimum solution. This refinement approach does have one implicit limitation, in that 

142 it cannot tell you directly the uncertainty of parameters within a model, only select 

143 the best set of parameters. In this study, we compare these best parameter sets, p. 
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144 for each misalignment and fook for variational trends in our parameters. This tests the 

145 sensitivity of each parameter to angular misalignment, provided the GA was successful 

146 at finding the global minima for the misaligned XRR data. The allowed parameter 

147 ranges in a model can sometimes influence the best solution, especially for intentionally 

148 misaligned data. Therefore, it is wise to use large ranges for parameter estimations in 

149 such a study. 

150 3. 3. Introduction of sample misalignment 

151 To simulate misalignment of the reflectometer zero angle, we shift the angle, 6'*, by a 

152 fixed amount, 6, 

154 for each angle in the initial, well-aligned dataset. We then generate simulated 

155 misaligned data which we then fit for each misalignment, 5 = [-0.025, -0.020, -0.015, 

156 -0.010, -0.005, 0, 0.005, 0.010, 0.015, 0.020, 0.025] in °. (Even though, for a typical 

157 XRR experiment covering an angular range of 3° from, a misalignment of 0.005° is just 

158 l/600th of the total range of the data, such a misalignment can have a substantial 

159 impact on parameter estimates, as we show below.) 

160 4. Results and Discussion 

161 The initial, well- aligned XRR data-set contained 620 {9j,I^) measurement pairs 

162 collected on a commercial XRR instrument at NMIJ using sample 1-08, covering a 

163 range of 0° to 3.1° in 9j in 0.005° steps. Each of the 11 data sets, one aligned, and 10 

164 misaligned, was fit using a GA refinement. Each GA refinement used a population of 

165 over 200 members with random parameters, p. Each GA was run for 5,000 generations 

166 and from several choices of initial parameter values in order to verify establishment of a 

167 true global minimum. Table [T] provides the allowed ranges for the 23 model parameters 

168 used in our refinements. Note that our prior knowledge of film thickness information 

169 was used to constrain multilayer thickness within a narrow window (1 nm for buried 

170 layers), and roughness to atomic feature scales (0.2 nm to 0.5 nm) assuming very-high 

171 quality interfaces with no interdiffusion; however, density was allowed to vary across a 

172 large range (0.5 to 2 times bulk density values). A second study was performed with a 

173 wider range of allowed roughness parameters (Table |2]) to address issues in refinement 

174 of surface layers for highly misaligned XRR data. We performed GA refinement with 

175 both NIST-developed software and a commercial package for comparison. Results from 

176 both software packages and both parameter ranges are shown. 

177 Figure [l] shows: a) the original data (aligned) GA refinement, and b) the 5 = 

178 +0.005° GA refinement. In both cases, the fit (line) to the data (points) are nearly 

179 perfect, with very little to indicate either misaligned or improperly modeled information. 

180 In comparison, figure [2] shows the results of: a) S = +0.025° GA refinement and b) 

181 6 = —0.025° GA refinement. In both cases, we see substantial deviations between 
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182 GA (line) and data (points). Several peaks are missed entirely in both refinements, 

183 and oscillation intensities vary substantially between model and data. This will have a 

184 significant impact on the density parameter estimates for these highly misaligned cases. 

185 Figure [3] shows the GOF results for the MSEiog GA refinements for NIST, 

186 commercial, and commercial wiht wide allowed roughness ranges. The better agreement 

187 between model and data from figure [T] is clearly seen by a minima in the range of 5 = 

188 to +0.005°. However, we do see a systematic shift between the GOF results from using 

189 Table [l] and Table |2] (wide). Wider roughness ranges accomodated more opportunities 

190 for surface roughness to exchange with surface layer thickness and provide lower minima 

191 GOFs for the negative offset angle, 5 < 0, cases. This illustrates the difficulty in using 

192 GOF as the mechanism for determining optimal alignment. Bias from allowed parameter 

193 ranges can directly affect GOF. 

194 Evaluation of top layer (s) thickness as a function of 5 led to uninterpretable results 

195 containing sharp discontinuities. Instead, we concentrated on layers deeper into the 

196 structure. Figure |4] shows thickness as a function of 5 for the bottom four buried 

197 GaAs and AlAs layers (see Table [l] for stack numbering). All layers show a surprisingly 

198 linear relationship suggesting the utility of using a simple linear extrapolation relating 

199 calculated thickness to misalignment. In Table [3] we show linear extrapolation results 

200 for all four layers. The table shows that sample 1-08 evaluated thickness and intercept 

201 show excellent agreement, indicating that the sample was extremely well-aligned. Note 

202 that the 5 = —0.025° was omitted in linear refinements due to non-linear behavior in 

203 some parameters. The slope term represents the change in thickness with change in 6. 

204 The negative slope indicates a decrease in thickness when data is misaligned in the +6^ 

205 direction. The values very close to 1 give us very high confidence in the linearity of 

206 the refinements. 

207 In figure [5| we see density as a function of 6 for the same four buried layers. 

208 Although all of the four layers show highly linear behavior, we see that the slopes vary 

209 depending on the layer under study. There is a trend for higher sensitivity in density 

210 determination with layers closer to the surface. In Table |4], we show both intercept 

211 and slope for linear fiting of the density variation data. In this case, we see several key 

212 features: First, the intercept densities from the model do not correspond to densities for 

213 either bulk AlAs or GaAs. For GaAs, the intercept density is higher, and for AlAs, the 

214 intercept is lower, than theoretical bulk values. This could either indicate a modeling 

215 bias, or it could tell us something about the density of GaAs and AlAs thin films grown 

216 by the molecular beam epitaxy process used, and therefore suggests an area for further 

217 study. Second, the slopes are both positive, and much larger (one order) than with 

218 thickness. From this we see that density increases with increasing shifts in 6^ and that 

219 density determination is highly dependent on, and sensitive to, sample alignment. This 

220 may explain some of the sample-to-sample and measurement-to-measurement variation 

221 in density often seen from XRR analysis. The R^ values for density are very close to 1, 

222 except in the case of layer 7 (AlAs next to GaAs substrate) which shows some variation. 

223 The linear nature of the thickness and density with respect to 6 suggests that 
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224 these linear extrapolations could provide us with a test for systematic errors of sample 

225 misalignment in commercial XRR instruments. Most XRR tools have an automated 

226 alignment procedure to validate the specular nature of a measurement prior to collection. 

227 Performing XRR scans with the GaAs/AlAs multilayer CRM through several mountings 

228 could establish the alignment repeatability for the instrument. In this scenario, some 

229 number of data sets (say 10) could be collected and analyzed with commercial GA 

230 software, and then the thickness information could be compared to the CRM values, to 

231 establish any systematic differences in thickness and anti-correlated changes in density, 

232 following our determined slope relationships. 

233 5. Conclusions 

234 This paper illustrates the proposed merits of using a Certified Reference Material, 

235 in this case a multilayer AlAs/GaAs CRM from NMIJ, in the assessment of sample 

236 alignment for the X-ray Reflectometry method. Sample alignment information can 

237 be better inferred through thickness (and possibly density) variation, than through 

238 more traditional GOF minimization. Over a reasonable range of sample misalignments 

239 (5±0.02°), both thickness and density vary in a linear way for the buried layers, allowing 

240 for direct extrapolation of misalignment effects. A combined strategy of assessing both 

241 thickness and density from buried layers of the multilayer structure may yield the 

242 best alignment assessment {U{5} < 0.005°). This combined approach needs further 

243 study, as density parameters are not certified to date. Note that this technique, by 

244 its nature, relies on the availability and validity of thickness and/or density certified 

245 parameters from a traceable reference material. By comparing systematic bias in these 

246 parameters, it is possible to assess potential misalignment errors caused by mounting 

247 and stage/sample alignment approaches on a commercial refiectometer. This is very 

248 important in the field, as sample alignment is often the dominant source of uncertainty 

249 in XRR measurement for high quality films. 

250 To date, these studies have only been developed using a single data set. Future 

251 work in this area will involve the automation of this alignment testing procedure to 

252 allow for repeatability testing of these results from multiple measurements and multiple 

253 instruments using the same specimen (CRM) to explore the limits of this approach. 
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Table 1. GA parameter ranges for multilayer refinement. 
Layer Material t/nm a/nm p/g cm~^ 

1 AI2O3 1.0-5.0 0.1-0.3 1.0-5.0 

2 GaAs 6.0-11.0 0.3-0.5 2.66-7.98 

3 AlAs 9.0-10.0 0.3-0.5 1.91-5.72 

4 GaAs 9.0-10.0 0.3-0.5 2.66-7.98 

5 AlAs 9.0-10.0 0.3-0.5 1.91-5.72 

6 GaAs 9.0-10.0 0.3-0.5 2.66-7.98 

7 AlAs 9.0-10.0 0.3-0.5 1.91-5.72 

8 GaAs - 0.3-0.5 5.316 



Table 2. GA with wide roughness ranges for multilayer refinement. 



Layer 


Material 


t/nm 


a/nm 


p/g cm ^ 


1 


AI2O3 


1.0-5.0 


0.1-2.0 


1.0-5.0 


2 


GaAs 


6.0-11.0 


0.1-1.0 


2.66-7.98 


3 


AlAs 


9.0-10.0 


0.1-1.0 


1.91-5.72 


4 


GaAs 


9.0-10.0 


0.1-1.0 


2.66-7.98 


5 


AlAs 


9.0-10.0 


0.1-1.0 


1.91-5.72 


6 


GaAs 


9.0-10.0 


0.1-1.0 


2.66-7.98 


7 


AlAs 


9.0-10.0 


0.1-1.0 


1.91-5.72 


8 


GaAs 




0.1-1.0 


5.316 





Table 3. Linear fit to thickness 


in layers 


from misaligned dati 


Layer 


Material 


Sample 


Intercept 


Slope 


coefficient of 






1-08 




m/ 


determination 






t/nm 


t/nm nm (°) ^ 




1 


AI2O3 


NA 








2 


GaAs 


NA 








3 


AlAs 


NA 








4 


GaAs 


9.281 


9.272 


-2.845 


0.996 


5 


AlAs 


9.468 


9.469 


-4.990 


0.999 


6 


GaAs 


9.269 


9.267 


-3.311 


0.997 


7 


AlAs 


9.458 


9.461 


-5.193 


0.994 


8 


GaAs 


NA 
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Table 4. Linear fit to density 


in layers 


from misaligned data. 


Layer 


Material 


Bulk 


Intercept 


blope 


coemcient or 










/ 

m/ 


determination 






p/nm 


p/nm 


nm (°) 


-1 T-)2 

^ R 


1 


AI2O3 


NA 








2 


GaAs 


NA 








3 


AlAs 


NA 








4 


GaAs 


5.316 


5.455 


47.24 


0.9997 


5 


AlAs 


3.81 


3.680 


38.05 


0.9995 


6 


GaAs 


5.316 


5.367 


37.53 


0.9998 


7 


AlAs 


3.81 


3.633 


19.01 


0.989 


8 


GaAs 


NA 










Figure 1. XRR measurement data and GA model refinement result for: a) aligned 
data and b) data intentionally shifted +0.005° 6 in w. Data is represented by points 
with a line representing GA refinement. 
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Figure 2. XRR measurement data and GA model refinement result for data 
intentionally shifted: a) +0.025° and b) -0.025° 6 in uj. Data is represented by points 
with a line representing GA refinement. Inset boxes are magnified views of two regions, 
showing deviations between refinement and data. 
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Figure 3. Goodness-of-fit for the GA refinement of the XRR multilayer model as a 
function of sample misalignment, S. NIST developed XRR code; is represented by x. 
A commercial GA XRR code is shown for comparison with + symbols. 
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Figure 4. Thickness determination of: a) GaAs layers, and b) AlAs layers, in 

multilayer stack, as a function of sample misalignment, 5. Note a decrease in slope 
for thickness as a function of sample misalignment. Dashed line represents the CRM 
thickness. 
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Figure 5. Density determination via genetic algorithm of: a) GaAs layers, and b) 

AIAs layers, in multilayer stack, as a function of sample misalignment, 5. Note the 
variation in the slope of density variance between layers in the stack. Layers 4 and 5 
show larger slopes (more sensitivity to misalignment). 



